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"^ : Abstract 

f~^ I The signal-to-noise ratios (SNRs) for quasi-circular binary black hole inspirals computed from 

CN ' restricted post-Newtonian waveforms are compared with those attained by more complete post- 

^~>' Newtonian signals, which are superpositions of amplitude-corrected harmonics of the orbital phase. 

It is shown that if one were to use the best available amplitude-corrected waveforms for detection 

templates, one should expect SNRs in actual searches to be significantly lower than those suggested 

'^ I by simulations based purely on restricted waveforms. 
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I. INTRODUCTION AND OVERVIEW 

Gravitational waves from the quasi-circular, adiabatic inspiral of compact binaries consist- 
ing of neutron stars and/or black holes are expected to be linear combinations of harmonics 
of the orbital phase, with the second harmonic as the dominant one. Such waveforms have 
been computed in the post-Newtonian (PN) approximation, where amplitudes and phases 
are expressed as expansions in the orbital velocity v (see [|| for a review and extensive ref- 
erences). For binaries without spin, the best PN waveforms currently available are of order 
v^ in amplitude and v"^ in phase [Q. In the usual notation this corresponds to 2.5PN and 
3.5PN order, respectively, and we will refer to such waveforms as being of (2.5, 3.5)PN order. 
Often the PN corrections to the amplitudes are discarded, in which case one ends up with a 
single harmonic at twice the orbital phase. This is known as the restricted post-Newtonian 
waveform. 

Searches for gravitational wave signals from inspiraling compact binaries are performed 
by matched filtering with a bank of templates. The same kinds of waveforms that go 
into template banks are also used as simulated signals injected into stretches of data to 
evaluate algorithms that search for real events and veto spurious ones. Several different 
types of waveforms are in use, all centered on the restricted post-Newtonian approximation. 
Apart from the straightforward restricted PN waveforms there are the Pade |^ and effective 
one-body ||^ waveforms, which result from resummation schemes designed to improve on 
the convergence of the PN phasing; these have no amplitude corrections either. Recently 
there has been much interest in the phenomenological templates proposed by Buonanno, 
Chen, Pan, and Vallisneri [P, which do have some simple corrections to the amplitude. 
Nevertheless, overall there has been a tendency to focus on phasing rather than amplitude, 
mainly because of the expectation that in matched filtering, the main issue is to know the 
phasing and number of cycles, or more precisely the number of useful cycles [0, of the signal 
in the detector's bandwidth. 

Here we investigate the effects of amplitude corrections at high PN order, or the absence 
thereof, on the reliability of templates and/or simulated signals. Our goal is to arrive at 
general statements concerning signal-to-noise ratios and detection rates. We will take both 
restricted and non-restricted waveforms to be the "standard" ones, so that the waveforms 
in one family are simple truncations of waveforms in the other, and the two families are 
parametrized analogously. Even so, when comparing waveforms belonging to different fam- 
ilies in the context of detection and parameter estimation, details regarding whether one 
compares faithfulness as above or effectualness are crucial [^, ^. In some of our considera- 
tions these are neglected; however, our main result will be largely independent of them. 

Given any two waveforms /ii, /i2, the signal-to-noise ratio (SNR) for the "detection" of 
/i2 using (the normalized counterpart of) hi as a template is given by 

r, 1 _ {.hi\h2) 

phAH = ^ ,, I, , ^ (1-1) 

y[hi\hi) 

where ( . | . ) is the usual inner product in terms of the noise power spectral density Sh{f) of 
the detector. With the convention x{f) = j_ x{t) exp(— 27rz/t) dt for Fourier transforms. 



The SNR for "detecting" a restricted PN waveform h^ using the same waveform for a 
template would be ph^lho], and this is the kind of SNR one encounters in past and cur- 
rent simulated searches. When using a restricted template bank in real data searches, one 
may expect the SNR for a genuine signal to be better approximated by Pho[h], with h an 
amplitude-corrected waveform. For simplicity, assume the parameter values of h and ho to 
be identical. As we will demonstrate, one then has 

pho[h] < pholho], (1-3) 

the difference being as large as 30% for astrophysically relevant sources. We stress that in 
the real detection problem one would maximize the SNR not just over extrinsic parameters 
but over intrinsic parameters as well. Thus the comparison of faithfulness in Eq. ( |1.3| ) is 
not adequate and one needs a comparison of effectualness @, |^ • To draw a firm conclusion 
further study is needed. 

However, the main focus will be on a different issue. The inequality ( [I.3[ ) refers to past 
and current searches. In future searches one might switch to amplitude-corrected templates. 
As we will show, one has ph[h] > Pho[h] with the two differing by only a small amount. 
Combining this with ( |1.3[ ), one arrives at 

Ph[h] < Pholho] (1.4) 

with the two differing by up to 25%. Hence, when using the best available amplitude- 
corrected waveforms for detection templates, one should expect SNRs in actual searches 
to be significantly lower than those suggested by simulations based purely on restricted 
waveforms. Since the associated detection rates are proportional to the cubes of the SNRs, 



they will differ by up to a factor of two. Note that unlike in the interpretation of (|1.3|) , no 
caveats are needed here (see the analytical treatment in Sec. |III B|) beyond the obvious fact 
that a genuine signal will still differ from a (2.5,3.5)PN waveform. 

In what follows we first outline the construction of amplitude-corrected waveforms in the 
stationary phase approximation. Next the above statements concerning SNRs are demon- 
strated numerically and explained analytically, after which we present conclusions. 

For concreteness we will work with the Initial LIGO design sensitivity from |^, but our 
results carry over to other initial detectors such as VIRGO. We set G = c = 1 unless stated 
otherwise. A waveform that is of PN order p in amplitude and q in phase will be called a 
(p, g)PN waveform. 

II. THE AMPLITUDE-CORRECTED WAVEFORMS 

The waveforms in the two polarizations take the general form 

h^,. = ^x {i/|% + xV2i7(V2) ^ ^^(1)^ ^ x'/'H^^;^^ + x2<i + x^/^i/^x^' } (2.1) 

where r is the distance to the binary, M its total mass, and r] the ratio of reduced mass to 
total mass. The post-Newtonian expansion parameter is defined as x = (27rMF)^'^ = v"^, 
with F{t) the instantaneous orbital frequency. The coefficients H^{^ , p = 0, ...,5, are 
linear combinations of various harmonics of the orbital phase with prefactors that depend 
on the inclination angle t of the angular momentum of the binary with respect to the line of 



sight as well as on 77; their explicit expressions can be found in 0. The measured signal also 
depends on the polarization angle if) and the position in the sky (6', 0) through the detector's 
beam pattern functions Fj^y. 

h{t) = F+h+{t) + F^h^{t). (2.2) 

For ground-based detectors, which are the ones we will be concerned with, it is reasonable to 
approximate -F+,x as being constant in time for the duration of the signal in the detector's 
bandwidth. The signal ( |2.2| ) is a linear combination of harmonics of the orbital phase \E'(t) 
with offsets (p(k,m/2)- We consider (2.5,3.5)PN waveforms, which contain seven harmonics: 

7 5 

^(^) =J2Y1 Ak,m/2){t) C0s(/C^(t) + iPik,m/2))- (2.3) 

fc=l m=0 

The index k runs over the harmonics while m/2 is PN order in amplitude. The harmonic at 
twice the orbital phase dominates, as it is the only one with a OPN amplitude contribution. 
In the restricted post-Newtonian approximation, no amplitude corrections are taken into 
account. In that case only one harmonic is present, namely the one at twice the orbital 
phase: 

ho{t) = A(2,o)(i) cos(2^(t) + (^(2,0)). (2.4) 

During the adiabatic part of the inspiral one has \d\nA(^k^s)/dt\ <^ kd'^/dt and 
\kd'^'^ /dt^\ <^ {kd'^ /dtY for k = 1,2, . . ., in which case one can use the stationary phase 
approximation (SPA) to the Fourier transform of (|2.3|) : 



hsPAif) = E 



fc=l 



ELo%W2)(Mi/))e~^-('=-/^) 



exp 

2.5 



2nft,-7r/4 + k^[^f 



^\/kF{t{lf)) 

(2.5) 

where [.J2.5 denotes consistent truncation to 2.5th post-Newtonian order (i.e., the "Newto- 
nian" prefactor /"^^^ is taken outside and the remaining expression is expanded in (27rM/)^/^ 
up to fifth power). A dot denotes derivation with respect to time and tc is the coalescence 
time. The function t{f) is defined implicitly through the instantaneous orbital frequency by 
F{t{f)) = f. The "frequency sweep" F in terms of F can be obtained from the expressions 
for energy and flux in [^, and the same goes for the phase ip{f)- 

Finally, the expressions of the respective harmonics are only valid below some upper cut- 
off frequency. The kth harmonic will be terminated at a frequency kfiso, where fiso is the 
orbital frequency corresponding to the last stable orbit. In the time domain this roughly 
corresponds to cutting off all of the harmonics at a time determined by F{t) = fiso- In the 
point mass limit one has Jlso = (6^/^27rM)~^, and for simplicity this is the expression we 
will adopt. To impose these restrictions on the harmonics, in practice we multiply the fcth 
harmonic by 9{kfLso ~ f)-, where 9{x) is the usual Heaviside step function. In the definition 
of the inner product ( . | . ), Eq. ( p..2|) , it is then natural to take fmax = Tfisoy the frequency 
reach of our highest harmonic. The lower cut-off frequency fmin is detector-dependent; for 
Initial LIGO it is set to 40 Hz. 

The SPA is just one method of approximating the Fourier transform; the main alternative 
is the fast Fourier transform (FFT). As shown in [0, ^ in the context of Initial LIGO, SNRs 
for restricted waveforms computed using the FFT and SPA differ by only a few percent for 



total masses up to ~ 30 Mq, and this is the mass range we will focus on. We also note that 
the conditions for the applicability of the SPA on the harmonics become more favorable with 
increasing k. 

III. SIGNAL-TO-NOISE RATIOS 

A. Numerical observations 

To set the stage, let us look at the SNRs in Initial LIGO for three different inspiral events 
at a distance of 20 Mpc: a binary neutron star (NS-NS), a neutron star and a black hole 
(NS-BH), and a binary black hole (BH-BH). The mass of the neutron star is taken to be 
1.4 Mq while for the black hole we choose 10 Mq. Table | displays the variation of pholho], 
Pho[h], and Ph[h] for (p, 3.5)PN waveforms h with increasing p, where Hq is the corresponding 
restricted waveform. The most interesting case is the NS-BH system. In going from p = 
to p = 0.5, both pho[h] and ph[h] increase, as has been noted for ph[h] by Sintes and Vecchio 
| rT| . However, at p = 1 they drop below p/j,j[/io], and a downward trend is seen as p continues 
further. At p = 2.5 one has Pho[h] < pho[ho], with pho[ho] being larger by 16.5% of Pho[h] - 
modeling signals as restricted waveforms and searching for them with restricted templates 
leads to an overestimation of SNR. Still at 2.5PN, Ph[h] > pho[h] but with ph[h] being larger 
by only 3.6% of Pho[h]] hence there is only a modest improvement to be had from the use of 
amplitude-corrected templates. Specifically, Ph[h] < Pho[ho] with a difference of 12.4%. The 
NS-NS and BH-BH systems exhibit similar trends. 
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Table I: Change in signal-to-noise ratios with increasing p in (p, 3.5)PN waveforms, for three dif- 
ferent systems at a distance of 20 Mpc, as seen in Initial LIGO. (Angles were chosen arbitrarily as 

6* = = 7r/6, 4^ = 7r/4, i = 7r/3.) 



In the Table we chose specific values for the angles in the problem. Once again consider 
the NS-BH system and let hhe a. (2.5, 3.5)PN waveform. Then for an optimally located and 
oriented source {9 = l = 0), Pho[ho] is larger than Pho[h] by 12.4%, but Ph[h] is larger than 
Pho[h] by only 0.14%. For a rather inconveniently located and oriented source (say, 6 = 7r/2, 
(f) = tp = 0, and L = it/2), the corresponding numbers are 19.2% and 5.1%, respectively. As 
these examples indicate, the relative difference between ph[h] and pho[ho] has only a weak 
angular dependence, varying between 12% and 14% for the system at hand. 

This behavior is quite generic. In the left panel of Fig. |I| we compare pho [h] , where h is a. 
(2.5,3.5)PN waveform, with ph^lho]. The fraction by which SNR is being overestimated by 
modeling signals as restricted waveforms is 5 = {pho[ho] —Pho[h])/Pho[h]- In a significant part 
of parameter space this fraction is between 20 and 30%. Indeed, on the line corresponding to 
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Figure 1: Plots of 5 (left panel) and 7 (right) as functions of rj and chirp mass Mc = M?]^'^. In 
both plots we have imposed that individual component masses be at least 1 Mq and that total 
mass M < 80 Mq; hence the sudden drop on the left in each panel. 

6 = 0.2, lighter and heavier component masses rrii, m2 are constrained to the ranges 1 Mq < 
mi < 12 Mq and 12 Mq < ma < 15 Mq, while for 6 = 0.3 one has IMq < mi < SM© and 
17 Mq < 7712 ^ 24 Mq. The implied range of binary systems is of clear astrophysical interest. 
As already suggested by the results of Table |, the overestimation is more pronounced for 
asymmetric systems and shows a strong dependence on mass. 

Next, let us look at the gain in SNR obtained by using amplitude-corrected waveforms, 
rather than restricted ones, as templates, with the signal being an amplitude-corrected 
waveform. The fractional gain is 7 = {ph[h] — Pho [h])/pho [h]- The right panel of Fig. |l| shows 
a plot of 7 as a function of chirp mass and 77. Clearly the use of non-restricted waveforms 
as templates would not give us much advantage over restricted ones; in most of parameter 
space the difference between the two is in the order of a few percent. Large gains (> 20%) 
do occur, but only for the astrophysically less interesting cases with large M^ and small 77, 
corresponding to very asymmetric binaries with high total mass (> 40 Mq). 

Roughly speaking, the left panel in Fig. |l] indicates by how much SNR has hitherto been 
overestimated due to inadequate signal modeling. The right panel shows that, even if in 
future searches one were to use amplitude-corrected waveforms throughout, this overesti- 
mation will not get compensated for. Indeed, putting together the information from both 
panels, we may conclude that 

Ph[h] < pho[ho] (3.1) 

with the two sides differing by up to 25% for sources with M < 30 Mq. 



B. Analytic considerations 

In this subsection, for brevity we focus on the inequality 

Pho[h] < Pholho], (3.2) 

but its analytic explanation will already provide some insight into the more interesting result 



( p.l|) ; a much more detailed analysis will be given in |jTO[. The above inequality is equivalent 



to {ho\h) < (ho\ho). Writing h as a sum of amplitude-corrected harmonics h^''\ the LHS 
may be approximated as 

7 

{ho\h) = Y,iho\h^'^)^{ho\h'^'^), (3.3) 

k=l 

because different harmonics tend to interfere destructively. (It can be checked numerically 
that the approximation is valid to within a few percent in the Initial LIGO band.) We stress 
that h^"^^ 7^ ho- Indeed, both are proportional to the second harmonic, but h^^' contains 
amplitude corrections while ho does not. 

The SPA for h^'^^ is of the form (leaving out the step function in frequency): 



h^'\f)=cr/' 



5^C(„/2)(27rM/)™/3 



m=0 



exp 



27r/te - 7r/4 + 2^ ( ^/ 



(3.4) 



where C is a real function of chirp mass and distance while the coefficients C(m/2) are complex 
functions of {9, 0, ip, t', v)- Setting C(^rn/2) = for m > 1 would lead to the restricted waveform 

Consider the form of the quantity Re[/iQ(/)/i*-^'*(/)] appearing in the integrand of (holh^'^^): 

5 

Re[h*{f)h^'\f)] = C'r"' Y. MChCim/2)]{2nMfr'\ (3.5) 

771=0 

This is to be compared with the expression Re[/io(/)/io(/)] = C^/^^/^|C(o)P in the integrand 
of (/lol^o). Now, it is not difficult to show that 

Re[qo)C(™/2)] < (3.6) 

for m = 1, . . . , 5, irrespective of parameter values. In particular, the functions Re[C/*gxC(i/2)] 
and Re[C,*QNC(3/2)] vanish identically while the others are negative definite. Consequently, 

MhlU)h^'\f)] < Re[K{f)ho{f)]. (3.7) 

Thus, one should expect {ho\h) < (/io|/io)) whence Pho[h] < pho[ho]- 

Next we revisit Table |. Let h he a {p, 3.5)PN waveform with < p < 2.5. If p = 0.5, 
only the first term in the RHS of ( p.5|) survives [|ll|. In that case {ho\h) > (/lol^o) due to 
the small contributions {ho\h^^^) and {ho\h^^^) to the LHS, which were neglected in ( p. 31) . 
However, as p is increased there is a notable downward trend in Pho[h] as a result of the 
inequalities (|3.6|) , except at p = L5 because of Re[C/Q^C(3/2)] = 0. 

The above analysis suggests that when using restricted templates it should be possible to 
construct adequate signal waveforms without having to include higher harmonics. Indeed, it 
should suffice to reinstate the amplitude corrections to the second harmonic only, in which 
case we arrive at the waveforms /i*-^-* . This is borne out in Table 0. 

Although slightly more involved, the inequality Ph[h] < Pho[ho] can also be proved ana- 
lytically by continuing along the lines above |10[. This direct proof not going via Pho[h] is 
more explicitly free of the issue of maximization over intrinsic parameters. 



NS-NS 
NS-BH 
BH-BH 



Ph,[h] Phoih^'n 



6.2462 6.2462 
11.5836 11.5731 
26.2166 26.2126 



Table II: Signal-to-noise ratios for detection of (2.5,3.5)PN signals h and the amplitude-corrected 
second harmonic h^'^' , using restricted templates, for three "canonical" systems in Initial LIGO at 
a distance of 20 Mpc. 

IV. CONCLUSIONS 

We have compared the performance of the best available amplitude-corrected post- 
Newtonian waveforms, namely the (2.5, 3.5)PN waveforms, with that of the restricted ones, 
both as templates and as simulated signals. Our results can be summarized as follows. 

(i) In simulated searches, templates and simulated signals alike have hitherto been mod- 
eled essentially as restricted waveforms. In real searches, a genuine signal will presumably be 
better approximated by an amplitude-corrected waveform, and we have seen that the associ- 
ated SNRs may be considerably smaller than those suggested by the simulations. However, 
issues concerning maximization over parameters were neglected; it would be of interest to 
perform a more in-depth study. 

(ii) More importantly, in the future one might use amplitude-corrected waveforms both 
for templates and simulated signals. In that case one may expect simulated SNRs to be 
much closer to the ones seen in actual searches conducted using the amplitude-corrected 
templates. We found that SNRs computed with amplitude-corrected templates and wave- 
forms are smaller than the ones with restricted templates and waveforms by up to 25% 
for astrophysically interesting sources, corresponding to a difference in detection rates by 
up to factor of two. This result is much more robust; except for the fact that a genuine 
signal will still differ from a (2.5, 3.5)PN waveform, it is independent of any other detection 
issues. It implies that if one were to use the best available amplitude-corrected waveforms 
for detection templates, one should expect SNRs in actual searches to be significantly lower 
than those suggested by simulations that involve only restricted waveforms. 

The effects seen here result largely from the amplitude corrections to the dominant har- 
monic. An interesting consequence is that by taking only these into account, one arrives at 
signal models which, given restricted templates, behave like the more complete waveforms 
despite the fact that they lack other harmonics. Overall, our results clearly underscore the 
importance of including amplitude corrections in templates and simulated signals. A more 
detailed treatment will be given in a forthcoming paper []TU|. 
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